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Abstract. Microscopic nuclear structure calculations have been performed within the 
framework of the unitary-model-operator approach. Ground-state and single-particle energies 
are calculated for nuclei around 14 C, 16 and 40 Ca with modern nucleon-nucleon interactions. 



1. Introduction 

Recently, ab initio nuclear structure calculations starting with a realistic nucleon-nucleon (NN) 
interaction have been possible beyond few-nucleon systems. In this kind of study, the microscopic 
derivation of an effective interaction from the realistic NN interaction allows one to perform the 
structure calculation in a restricted model space. One of the successful methods of the structure 
calculation in this direction would be the no-core shell model (NCSM) ^|2j. In the NCSM, the 
microscopic effective interaction (Hamiltonian) that takes account of the short-range correlation 
of the original NN interaction can be derived through a unitary transformation of the original 
Hamiltonian. Although the shell-model diagonalization is done in a large model space so that 
the final results do not depend on the model-space size, such a model space is still much smaller 
than the huge Hilbert space of the original Hamiltonian. The microscopic understanding of 
nuclear structure has been growing through the NCSM as well as the Green's function Monte 
Carlo (GFMC) E] in which the original NN interaction and a three-nucleon force can be 
directly used. However, although the GFMC and the NCSM are powerful method to describe 
nuclear structure, the application of these methods to the nuclear structure calculation may be 
limited to light nuclei up to A ~ 12 due to the present computer power. 

If one wishes to describe heavier nuclei, one needs to have another method. One of the 
promising methods may be the unitary-model-operator approach (UMOA) [310 [7j. The UMOA 
can be regarded as one of the coupled-cluster methods (CCM) [SJ E3 IH1 E] of Hermitian 
type. As for recent developments of the CCM in nuclear theory, one may refer to Refs. |13U14j . 
In the UMOA, an energy independent and Hermitian effective interaction is derived through a 
unitary transformation of the original Hamiltonian [151 116) which is essentially the same as the 
unitary transformation used in the NCSM. By doing the unitary transformation in a two-step 
procedure, the structure calculation can be performed beyond p-shell nuclei. So far, we have 
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performed structure calculations for not only stable nuclei around O but also neutron- 
rich oxygen isotopes ^Tj\ and A hypernuclei 0. In the following sections, we shall outline the 
calculation method of the UMOA and present recent results for nuclei around 14 C, 16 O and 
40 Ca. 



2. Method of calculation 

In the UMOA, the Hamiltonian H to be considered is obtained through a unitary transformation 
of the original Hamiltonian H as H = e~ s He s . The exponent S is an anti-Hermitian operator 
and written as S = arctanh(a; — u>) [SJ with a mapping operator ui = QluP under the 
restrictive conditions PSP = QSQ = 0, where P and Q are the usual projection operators and 
have the properties as P + Q = 1, P 2 = P, Q 2 = Q and PQ = QP = 0. Thus, the operator 
ui satisfies the relation lo 2 = uj^ = 0. We should note here that the unitary-transformation 
operator U = e s is also given by a block form concerning the P and Q spaces as 

TT-l P{l + ^Ujy l / 2 P -P W t(l +ww t)-V2Q\ 

{ Quil + u^-Wp Qil + uu^-WQ ) [ ) 

which agrees with the unitary transformation by Okubo 15 j. By applying the above unitary 
transformation to a two-body subsystem of the original Hamiltonian, the two-body effective 
interaction vi 2 of Hermitian type is given by 

v 12 = U-\h + v 12 )U -h , (2) 

where ho is the one-body part and v± 2 is the bare two-body interaction. The effective 
interaction in the P and, if necessary, Q spaces is determined by solving the decoupling equation 
Qv\ 2 P = PvyiQ = 0. The actual method of calculating the matrix elements of U and v\i with 
the harmonic-oscillator (h.o.) basis states in the neutron-proton (np) formalism may be found 
in Ref. 0. 



2.1. Two-step calculation for the effective interaction 

In order to make the structure calculation in an inexpensive way, we perform the unitary 
transformation twice as follows. First, we derive the two-body effective interaction in a large 
model space to take into account the short-range correlation of the original NN interaction. The 
large model space consisting of two-body states is specified by a boundary number p\ which is 
given with the sets of h.o. quantum numbers {n a ,l a } and {nb,h} of the two-body states by 
pi = 2n a + l a + 2nb + lb- The value of p\ is taken as large as possible so that the calculated results 
do not depend on this value. If we diagonalize the transformed Hamiltonian with the many-body 
shell-model basis states using the effective interaction in the large model space, that leads to the 
NCSM. However, the calculations for heavier nuclei such as 16 and 40 O may not be practical 
in this manner because the present computer cannot diagonalize the huge matrix elements of 
the Hamiltonian which is large enough to obtain the converged results. If we intend to obtain 
only the energies of the ground states of closed-shell nuclei and the single-particle (-hole) states 
in its neighboring nuclei, it would be convenient to perform the unitary-transformation again 
as follows. We define a small model space Pnp and its complement Qnp by separating the large 
model space in the previous procedure as shown in Fig.^ The symbols p n and p p in Fig.^stand 
for the uppermost occupied states of the neutron and proton, respectively. We here only show 

(2) (2) 

the np channel. It should be noted that the Pnp and Q n p spaces are considered on an equal 
footing in this second-step calculation when we solve the decoupling equation for the effective 

(2) (2) (2) (2) 

interaction v\ 2 in this step as Qnp v\ 2 Pnp = using the effective interaction determined in the 
first-step calculation as an input. Namely, we derive the effective interaction again in the large 
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Figure 1. The spaces consisting of two- 
body states for the np channel considered in 
the second-step calculation for the effective 
interaction. 



model space. However, by taking the model space Pnp and its complement Qnp as shown in 
Fig. ^ the resultant effective interaction in this second step has no vertices which induce 
two-particle two-hole (2p2h) excitation. This is analogous to the Hartree-Fock (HF) condition 
which means that an original Hamiltonian is transformed so that the matrix elements for lplh 
excitation reduce to zero. Although the vertices of the one-body non-diagonal matrix elements 
remain in determining the effective interaction in the UMOA, these non-diagonal matrix elements 
are diagonalized at the end of the structure calculation. We notice that the two-body effective 
interaction in the first- and second-step calculations is determined self-consistently with the one- 
body potential for both the particle and hole states. Although this procedure is not necessarily 
needed in determining the effective interaction, such an effective interaction may be optimized 
for a restricted model space so as to obtain a good unperturbed energy. 

We here show an example of the effect of renormalizing 2p2h excitation in the second-step 
calculation for the CD-Bonn potential |18| . The unperturbed ground-state energies of 16 for 
pi = 12 and htt = 15 MeV using the effective interactions determined in the first- and second- 
step calculations are —38.40 and —104.12 MeV, respectively. This means that a large amount 
of the 2p2h effect is renormalized into the effective interaction in the small model space in 
the second-step calculation. In order to obtain the final result including the rest correlation 
effects, we diagonalize the transformed Hamiltonian with the shell-model basis states, taking 
into account lplh excitation from the unperturbed ground state. As for the closed-shell nucleus 
plus one-particle (one-hole) system, the shell-model basis states are composed of the lp and 2plh 
states (lh and lp2h states). These correlation energies are added to the unperturbed ground- 
state energies, and then we obtain the final results of the ground-state energies of the closed-shell 
nucleus and the closed-shell nucleus plus one-particle (one-hole) system. We note here that there 
remain some correction terms to be evaluated, such as the three-body cluster terms. Although 
the three-body cluster effect is essentially small, that may have a significant effect for some 
particular cases. Actually, we have found that the evaluation of the three-body cluster effect 
plays an important role to describe shell structure in neutron-rich oxygen isotopes However, 
generally speaking, the magnitude of the three-body cluster effect is much smaller than the two- 
body cluster effect, and thus the cluster expansion of the transformed Hamiltonian is justified 
from the viewpoint of the perturbative expansion. In the next section, we show some of the 
recent results which do not include the three-body cluster effect. The results including the 
three-body cluster effect will be reported elsewhere in the near future. 



Table 1. The calculated ground-state energies with the lplh effect E gmSm and the binding 
energies per nucleon BE/A of 16 0. All energies are in MeV. 
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3. Results and discussion 

First of all, we show calculated ground-state energies of 16 with the lplh effect with modern 
NN interactions in Tabled The final results for the Nijm-93, Nijm-I |19j . the CD-Bonn ^H] 
and the chiral N 3 LO [2U] potentials are tabulated together with the experimental values. The 
binding energies per nucleon are also shown. The Coulomb force is used commonly. The results 
for the Nijm 93 and the CD Bonn are the least and most attractive, respectively, of the four 
potentials. This tendency can also be observed in the Faddeev-Yakubovsky calculations for 4 He 
by Nogga et al [U- It is seen that the calculated ground-state energies are less bound than 
the experimental value. In the present results, higher-order correlations such as the three-body 
cluster terms have not been evaluated. If we include the effect of the three-body cluster terms, 
the calculated results become more attractive by a few MeV. For example, in the case of the 
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Figure 2. The calculated single-particle and single-hole energy levels in O and O, 
respectively, with modern NN interactions relative to the ground state of 16 0. 
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Figure 3. The p\ and Tiil dependences of 
the calculated ground-state energies of 40 Ca 
for the CD-Bonn potential. 




Figure 4. The calculated single-particle 
energies in 15 C for the CD-Bonn potential. 



CD-Bonn potential, the three-body cluster effect is about —4 MeV. However, even if we add this 
value to the ground-state energy in Tabled the result is less attractive than the experimental 
value. In the present calculation, a genuine three-body force is not included. The inclusion of 
the real three-body force and the higher-order many-body correlations would compensate for 
the discrepancies between the experimental and calculated values. Such a study remains as an 
important task for a deeper understanding of nuclear ground-state properties in the UMOA. 

In Fig. |21 calculated single-particle levels in 17 and single-hole levels in 15 relative to the 
ground state of 16 are illustrated. The values in Fig. |2 are the spin-orbit splitting energies for 
hole states in 15 and particle states in 17 0. We see that the calculated spin-orbit splittings for 
the hole states in 15 are smaller than the experimental value. As for the spin-orbit splittings 
for the particle states 17 0, as opposed to the hole-state case, the energies are larger than the 
experimental value. This may be due to an insufficient treatment for the 3/2 + unbound state 
in the present calculation. We have used only the h.o. states as the basis states. However, 
we may say that the calculated spin-orbit splittings for the hole and particle states in nuclei 
around 16 are, on the whole, not very different from the experimental values though the results 
somewhat depend on the interactions employed. The magnitudes of remaining discrepancies 
may be reduced if we include the genuine three-body force in the calculation and evaluate the 
higher-order cluster terms. 

We are now trying to do calculations for heavier systems. Here we show some of the 
preliminary results for 40 Ca. Figure 01 shows the p\ and hQ dependences of the ground-state 
energy with the lplh effect for the CD-Bonn potential. It is seen that a fairly convergent result 
is obtained at p\ = 16 and TiVt = 14 MeV. Though we are now calculating for p\ = 18, the 
difference of the results between p\ = 16 and p\ = 18 would be a few MeV at most. In the case 
of 16 0, the convergent result can be obtained at p\ = 14. Since we are now calculating heavier 



systems, we need a larger model space to obtain the convergent result. 

In Fig. 0J we show calculated single-particle energies in a neutron-rich nucleus 15 C together 
with the experimental values. The result for "unpt." denotes the unperturbed single-particle 
energy which is defined as the h.o. kinetic energy plus the self-consistent one-body potential 
determined in the second-step calculation. In 15 C, the ordering of the experimental single- 
particle l/2 + and 5/2 + states are opposite to the case of O. However, our result for "unpt." does 
not reproduce this tendency, and the single-particle levels are rather repulsive to the experimental 
values. However, if we see the results with the 2plh effect by the diagonalization as shown in 
the middle of Fig. 21 the results become more attractive and the two levels are reversed, and 
then a good agreement with the experimental value is obtained. We should remark, however, 
that the three-body cluster terms remain to be evaluated. As has been reported for neutron- 
rich oxygen isotopes in Ref. |17j . the three-body cluster has a non-negligible effect on particle 
states of loosely-bound neutron-rich systems. Therefore, for the complete study, we have to 
evaluate the three-body cluster terms and also investigate the effect of the genuine three-body 
force though the present two-body cluster approximation should be a good treatment. 

As has been shown before, the UMOA is a useful many-body theory to microscopically 
describe nuclear structure near closed-shell nuclei. The ground-state energy and the single- 
particle (-hole) energy can be calculated systematically for not only N ~ Z nuclei but also 
neutron-proton asymmetric systems beyond p-shell nuclei. Forthcoming new facilities for 
accelerating RI beams in the world will reveal new and exciting phenomena such as the change 
of shell structure in asymmetric nuclei. The UMOA has a possibility to develop new structures 
from a microscopic point of view in investigating the exotic nuclear systems. 
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